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ABSTRACT 

Aims. We study the individual evolution histories of three nearby low-mass edge-on galaxies (IC 5052, NGC 4244, and NGC 5023). 
Methods. Using resolved stellar populations, we constructed star count density maps for populations of different ages and analyzed 
the change of structural parameters with stellar age within each galaxy. 

Results. We do not detect a separate thick disk in any of the three galaxies, even though our observations cover a wider range in 
equivalent surface brightness than any integrated light study. While scale heights increase with age, each population can be well 
described by a single disk. Two of the galaxies contain a very weak additional component, which we identify as the faint halo. The 
mass of these faint halos is lower than 1% of the mass of the disk. 

The three galaxies show low vertical heating rates, which are much lower than the heating rate of the Milky Way. This indicates that 
heating agents, such as giant molecular clouds and spiral structure, are weak in low-mass galaxies. 

All populations in the three galaxies exhibit no or only little flaring. While this finding is consistent with previous integrated light 
studies, it poses strong constraints on galaxy simulations, where strong flaring is often found as a result of interactions or radial 
migration. 

Key words, galaxies: spiral, galaxies: individual: IC 5052, NGC 4244, NGC 5023, galaxies: evolution, galaxies: stellar content, galax¬ 
ies: structure 


1. Introduction 

Studies of galaxy evolution have always had to solve the problem 
of the very long timescales involved: Galaxies only change very 
slowly over millions and billions of years. The study of galaxy 
evolution must therefore rely on indirect methods and measure¬ 
ments. 

Fortunately, the finite speed of light opens a way to look into 
the past. By observing galaxies far away, we automatically look 
at galaxies in the past. This allows us to study the evolution of 
galaxy properties from the early Universe to the present time. 
This leads to a good knowledge of the evolution of the overall 
galaxy populations, but it does not help with the study of the 
evolution of individual galaxies. The question what the prede¬ 
cessor of the Milky Way, or of any current galaxy in general, has 
looked like cannot be answered through such studies. 

A different approach to studying galaxy evolution is the de¬ 
tailed analysis of current galaxies in order to distinguish their 
individual histories. The detailed knowledge of the kinematics, 
ages, and chemistry of all stellar populations of a galaxy today 
would allow inference on its formation and evolution. This ap¬ 
proach is often called galactic archeology or near-field cosmol¬ 
ogy. But even without the full knowledge of the kinematical and 
chemical distribution of the stars, galaxies contain much infor¬ 
mation about their past. Integrated colors (and their gradients) 
contain information about the ages and metallicities of the un¬ 


derlying stellar distribution, integrated spectra give valuable in¬ 
formation about the ages, kinematics, and metallicities and even 
the temporal evolution of them (e.g. Qcvirk et al.| 2008), and 
deep color-magnitude diagrams allow determining the star for¬ 
mation history (SFR) and metal enrichment function. 

Most processes in the evolution of galaxies also leave their 
signature on the structure of galaxies; an active merger history, 
for instance, will lead to a strong bulge, while the quiescent ac¬ 
cretion of gas forms a thin galactic disk. In this paper we aim to 
better understand galaxy evolution by analyzing the vertical and 
radial structure of stellar populations in nearby galaxy disks. 


1.1. Vertical structure of disks 

From simple modeling of the vertical structure of galactic disks, 
it is expected that it follows a sech^ density profile when an 
isothermal disk is assumed ( |van der Kruit & Searle|1981a| , but 
observationally, a centrally more peaked profile such as a sech 
or exponential profi lj^was found t o be more appropriate for the 
stellar distribution ( [van der Kruit||198^ |de Grijs et al.||1997] ). 
This deviation from the isothermal model can be explained by a 
mixture of stellar populations with different velocity dispersions, 

^ Often the more general sech^/”(j7z/zo) is used, which also includes 
an exponential forn ^ oo. These functions differ only in the midplane, 
they are all asymptotically exponential for large z. 
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as they are observed in the Milky Way, where the velocity dis¬ 
persion increases with stellar age (|Wielen| 19771 |Carlberg et aL| 
IT9851 ). 

Observations show that a single-disk profile is often not suf¬ 
ficient to describe the disk, but that a second component with a 
larger scale height is necessary. These so-called thick disks wer e 
first discovered in SO galaxies ( |Burstein|1979t |Tsikoudi|1979|), 


then in many other galaxies and in our own Milky Way ( [Gilmore 
& Reid||1983|). Later they were found to be ubiquitous ([Pohlen 


et al.|2004||Yoachim & Dalcanton|2006l|Comer6n et al.|2011a| ). 

Thick disks were usually seen as a distinct component. Many 
properties of the thick disk could be determined in the Milky 
Way: It is kinematically hot (e.g., [Nissen 1995[ [Girard et al 


2006), old, metal poor, and alpha enhanced (e.g., jFuhrmann 


2008| ). Recent ly, the picture of a clea r distinct thick disk has been 
challenged by Bovy et al. ( |2012a|b| ). They proposed the disk to 


be a superposition of many “mono-abundance population disks”, 
which can be each described by a single-disk model. 

In many studies it was found that thick disks do not only have 
larger scale heights, but also larger scale lengths ([Robin et al. 


[1996t[^achim & Dalcanton|2006t[Juric et al.[20'08l and others). 


but thi s ha s also been question ed for the MW by [Bensby et al. 
( 201 1[) and[Bovy et al. ([2012b[) as we ll as by recent simulations 


2009 


Loebman et al 


dial migr 

MM 


^ b ut this effect was s hown to be too weak to form the observed thick 
disks ( Villumsen|1985[ ) 


We note that various more recent simul ations ( [Minchev et al.|201^ 


[Martig et al.[2014b[[Vera-Ciro et al.|20T4] l have cast doubt on the heat 
ing effect of radial migration. 


an exponential light profile is a still-unsolved problem. The two 
prevailing ideas are that (1) the exponential profile refiects the 
initial angular mome ntum distribution of the gas cloud ( [Freeman 


1970 ; [L^son[[ 197 6) and (2) the viscosity of the gas leads to 


redistribution of angular momentum that results in an exponen 
tial pro file ([Lin & Pringle[[1987] ). More recently, Elmegreen & 


Struck[ ( [2013[ ) suggested that stellar scattering off of transient 


mass clumps in the disk naturally results in an exponential pro¬ 
file. 

A widely accepted idea is that disks form from inside out: 
at first, gas with low angular momentum gathers in the center 
and forms stars, and later gas with higher angular momentum 
gathers around this and forms more stars. This behavior is seen 
in many different models (e.g., Larson][1976[ White & Frenk 
[Burkert et al.[[1992[ [Mo et al.[[1998t [Naab & Qstriker 


1991 


2006 


Brook et al.|2006[) ^d leads to negative age and metallic- 


ity gr^ents (e.g.,[Matteucci & Francois|1989t[Chiappini et al. 

[1997t [Boissier & Prantzos[[1999t [Prantzos & Boissier[[2000[ ). 
Observational evidence for the inside-out formation of disks can 


( Stinson et al.|2013t[Bird et al.|2013|[Minchev et al.|2014[ ). 

The possible formation processes of the thick disk are still 
unclear. Thick disks may 

- form as thick disks, for example, in very massive star for- 
mation aggregates ([Kroupa|2002 ) or turbulent clumpy disks 
( Bournaud et al.[[2009[) suc h as would appear after gas-rich 
mergers ( [Brook et al.|20()4| ); 

- be created from a pre-existing thin disk that is thickened 
through internal heating or scattering by giant molecular 
clouds or spiral arm^ 

- be created from a pre-existing disk by redistributing stars 
through outward radial m igration ( [Schonrich & Binney[ 


Mac Arthur et al.|2009 

Williams et al.|2009[), in metallicity gra- 

dients (Mac Arthur et a 

[.|2009[), color gradients ([de Jong|1996[), 

and in the change in ga 

axy sizes with redshift ([Patel et al.|2013l 

[van der Wei et al.|2014 



The interpretation of radial gradients has to take into account 
that stars can change their radial position in the galaxy with time. 
It has long been known that scattering effects can blur and heat 
galactic disks. But it has been found only relatively recently that 
stars can also change their radial position and still keep the circu¬ 
lar nature of their orbits ( [Sell wood & Binney|2002[ ). This radial 
migration process dramaticall y changes the chemical evolu tion 
and the metallicity gradients ( [Schonrich & Binn^[20Q9[) and 
leads to a radial extens ion of disks ( [Roskar et al.|2008|[Sanchez-[ 
Blazquez et al.[ 2009[ ). It might also cause the formation of a 


- result from external heating of an earlier disk by minor merg¬ 
er^ Qum^£rar[|T99^ or dark matter halo bombardment 
( [Kazantzidis et al.|2008[ ); 

- be made of accreted material from satellite galaxies ([Staffer 
[T98^[Abadi et al.|2QQ3l . 

Obviously, a combination of these processes might also play a 
role, with different ratios in different types of galaxies. 

1.2. Radial structure of disks 

Disks are thin, rotationally supported structures. They show a 
radial light profile that is (close to) exponential ( [Patterson|1940[ 
[de Vaucouleurs|195^ [Freeman] 1970[). However, some disks are 
truncated in the outer parts ( van der Kruit[[1979[ ), which means 
that outside the break radius (typically at a radius of 1.5-6.0 scale 
lengths) a different exponential is needed to describe the surface 
brightness profile. In most cases (60%) the outer exponential is 
steeper than than the inner one, but for 30% it becomes shal¬ 
lower. Only a minority of 10% does not show a break in the 
radial light profile at all ( [Pohlen & Trujillo|2006[ ). 

Disks form through the dissipational collapse of a cooling 
gas cloud that forms stars after collapsing. How this leads to 


thick-disk component ( [Loebman et al.|20lT] ). 

While the physical reason for the exponential disk profiles 
is still unclear, the reasons for a break in the radial profiles are 
even less well understood. Th e breaks occur at a similar sur¬ 
face brightness in all galaxies ( [Kregel et al.[[2004[ ), at the same 
radius for all heights above the plane, and for all populations 
( Pohlen et al.|20Q7t[de Jong et al.|2QQ7[). They can be observed 


even at redshifts of z ^ 1 ( [Perez|2()04[ [Trujillo & Pohlen|2005[ ). 

Truncated disks have a color minimum at the break radius, while 
antitruncated, that is, disks whose scale length is greater be¬ 
yond the break, and unbroken disks only show a fiattening of 
their color profile at large radii, both in the local Universe and 
at higher redshifts (z 1) ([Bakos et al.[[2008l [Azzollini et al. 

|2008a|b[ ). -- 

Th e nature of breaks is also still discussed. IBakos et al.l 
( [2008[ ) claimed that a break is only due to a change in stellar 
populations and that the mass profile is unbroken, but this is in 
contrast to results from the GHOSTS survey ( de Jong et al. 2007] 
[Radbum-Smith et al.|2012[ ), which show breaks in star counts of 
different populations. The break in simulations is also connected 
to a steepening of the stellar mass profile ( [Roskar et al. [[20081 
Sanchez-Blazquez et al.|2009 ). These simulations predict a min¬ 


imum of the age distribution at the break radius and a smooth 
metallicity profile. 

The numerous break formation models can be roughly di¬ 
vided into two groups. The first connects the break to a break in 
star formati on, either due to a limited gas di stribution ([van der 


Kruit[[1987[) or to a star formation threshold ( Fall & Efstathiou 


1980[ ^nnicutt[[1989t [Dopita & Ryd^[1994[ Schaye 2004[ 


Elmegreen & Hunter|20b6 ). The second group contains models 
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that dynamically redistribute stars after their formatio n, either as 
a result of secular angular momentum redistribution ([Debattista 


et al.|2006} |Foyle et al.|2008[|Minchev et al.|2012b| ) or of tidal 

interactions ( |Gnedin|20()3t|Kazantzidis et al.|2008| ). 


We here examine the structure of different stellar populations 
with distinct ages to study the temporal evolution of galaxy 
disks. In Sect. we explain the data and methods, we present 
our results in Sect. [3l discuss them in Sect.|^ and conclude with 
a summary in Sect.^ 


2. Data and methods 

2.1. GHOSTS Survey 

GHOSTS is an extensive multicycle HST survey that images the 
resolved stellar populations in the halos and outer disks of 17 
nearby disk galaxies. The galaxies in the GHOSTS sample span 
a wide range of morphologies, from Sab to Sd, and masses (with 
Vrot ranging from 80km/s to 230km/s). A detailed description 
of the survey is given in |Radburn-Smith et al. ( 2011 ) and the 
data are publicly available at the Mikulski Archive for Space 
Telescopes (MAStQ 

The survey was performed using the cameras ACS and 
WFC3 and filters F814W and F606W. These filters were cho¬ 
sen because they have a high throughput and because red giants, 
which are supposed to be the majority of (bright) stars in the as¬ 
sumed old halos, have the highest flux in this wavelength range. 

The observations were designed so that red giant branch stars 
could be well resolved; the observations typically reach at least 
a signal-to-noise ratio (S/N) of 10 at 1 mag below the tip of the 
red giant branch. For the nearer galaxies (D < 5 Mpc) this re¬ 
quirement was obtained within the limits for HST SNAP pro¬ 
grams, while for farther galaxies GO programs were run. The 
dedicated GHOSTS observations were complemented with all 
archived observations fulfilling the same requirements. The cur¬ 
rent GHOSTS database contains data from cycles 12 to 21. 


the culls were requested to maximize the recovery fraction of 
artificial stars that are placed into the same images. We call the 
culls we chose with those two requirements sparse field culls be¬ 
cause they are optimized for fields with low star count numbers. 

A second set of culls was optimized for a high recovery frac¬ 
tion of artificial stars in very crowded regions. This is called the 
crowded field cull. Compared to the sparse field cull, it has about 
60% more contaminants, but at the same time, the recovery frac¬ 
tion of stars in crowded regions is more than twice as high as 
for the sparse field culls. Since the typical number of contami¬ 
nants is very small (about a few dozen per ACS field) compared 
to the number of stars in crowded regions (a few thousand per 
ACS field), the larger number of real detections outweighs the 
increased number of contaminants. For more details on deter¬ 
mining optimal culls we refer to the appendix of |Radbum-Smith 
et al.| ( MTT] ). In the following we always use the crowded field 
culls, if not stated otherwise. 

To further in crease the reliability of o ur star catalogs, we 
used SExtractor ( Bertin & Arnouts||1996 ) to create a segmen¬ 
tation mask that masks out all extended sources such as back¬ 
ground galaxies or nearby bright stars. 

2.1.2. Sample of edge-on galaxies 

The GHOSTS sam^e contains seven edge-on galaxies that have 
data on their disk^ which we included in our analysis (see 
Table[^. In this first paper we examine the three low-mass galax¬ 
ies (75km/s < Vrot < lOOkm/s). The three massive galaxies 
(Vrot > 200km/s) will be covered in a subsequent paper. The 
intermediate-mass galaxy NGC 4631 is analyzed in an additional 
paper because the strong effects of the interaction with its neigh¬ 
bor NGC 4627 require a different treatment. 


2.2. Age selection in the CMD 

Color-magnitude diagrams (CMD) allow us to separate the stars 
into populations of different ages. The CMDs of the GHOSTS 
galaxies show clearly distinguishable structures that belong to 


2.1.1. GHOSTS data reduction 

All detai ls of the observations and th e data reduction can be 
found in |Radburn-Smith et ak] ( |2Q11[ for the ACS data) and 
Monachesi et al. ^ 2015 for the WFC3 data). Here, we only sum¬ 
marize the main points of the data reduction. 

After basic data reduction (bias subtraction, fiat-fielding, 
cosmic ray and bad pixel identification, drizzling), the program 
DOLPHOT, which is a modified version of HSTphot ( |Dolphin| 
2000| ), is used to identify and photometrically measure the stars. 
We performed PSF photometry using the tiny tim point-spread 
functions (PSF). Together with the positions and magnitudes of 
the stars, DOLPHOT reports many photometric parameters for 
diagnostic purposes, such as sharpness, roundness, crowding, 
and the S/N. These parameters can be used to distinguish be¬ 
tween stars and other objects. 

We defined a set of culls on the diagnostic parameters that 
ensured that a maximal number of stars and a minimal number 
of contaminants was detected. To estimate the number of con¬ 
taminants that are expected for a given set of culls, observations 
from the Hubble archive were chosen that are aimed at high- 
redshift objects at high galactic latitudes. These observations are 
expected to be free of any detectable stars, and thus the culls 
should minimize the detections in these fields. At the same time. 


different stellar populations (see the CMDs in Figs. to B3 


the appendix). These populations are the upper main sequence, 
the helium-burning branches, the asymptotic giant branch, and 
the red giant branch. 

To define age-separated stellar populations, we defined five 
CMD bins (see Fig.[^ left). These CMD bins were designed to 
cover different ages with as little overlap in age as possible (see 
Fig.[^ right). To define the bins and determine their age distribu- 
tions, we used a synthetic CMD created with MA TCH ([Dolphin 


1997 


2010 


2002|) and using th e Padova isochrone set ( [Girardi et al. 
Marigo et al.|2008| ). stant star formation rate (from 
log(t)=10.15 to log(t)=6.6) and a fiat metallicity distribution 
function (between [Z]=-2.2 and [Z]=0.2) were used. 

The populations defined by these CMD bins are as follows. 
Main sequence (MS, blue): contains massive core hydrogen¬ 
burning stars that are mainly younger than 40 Myr, but stars 
up to 300 Myr can contribute to this population. 

Upper helium-burning branches (upHeB, cyan): contains mas¬ 
sive core helium-burning stars with an age of mainly be¬ 
tween 40 Myr and 150 Myr, with a weaker contribution from 
stars between 25 Myr and 300 Myr. 


http://archive.stsci.edu/pub/hlsp/ghosts/ 


^ An eighth galaxy, NGC 5907, turned out to be farther away (D= 
16.8 Mpc) than previously measured. Therefore its CMD is not deep 
enough for our analysis, and it was excluded from further observations 
due to the long exposure times that would have been needed. 
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Table 1. Details of the galaxy sample. 


name 

RA(J2000.0)' 

(°) 

DEC(J2000.0)‘ 

(°) 

Zo' 

(") 

PA^ 

D 

Incl.^ " 

(°) 

morph. 

V ^ 

^ max 

(km/s) 

(m-M)o^ 

(mag) 

disF 

(Mpc) 

IC5052 

313.006792 

-69.19331 

13.35 

-38.0 

90.0 

7.1 

79.8 

28.76 

5.6 

NGC5023 

198.052500 

44.041222 

9.29 

27.9 

90.0 

6.0 

80.3 

29.06 

6.5 

NGC4244 

184.373583 

37.807111 

22.14 

42.2 

88.0 

6.1 

89.1 

28.21 

4.4 

NGC4631 

190.533375 

32.541500 

13.74 

82.6 

85.0 

6.5 

138.9 

29.34 

7.4 

NGC 891 

35.639224 

42.349146 

11.86 

22.8 

88.0 

3.0 

212.2 

29.80 

9.1 

NGC7814 

0.812042 

16.145417 

- 

134.4 

70.6^ 

2.0 

230.9 

30.80 

14.4 

NGC 4565 

189.086584 

25.987674 

13.67 

-44.8 

90.0 

3.2 

244.9 

30.38 

11.9 


Seth et al. 


( |2005a i;^^^ po sition angle, inclination, mor phology index 


j201lt . 


Notes. Right ascension, declination, an d scale heights froni K-Band fits by 

and rotational velocity from HyperLEDA ( [Paturel et al. 20031; distance modulus and distance from Radburn-Smith et al.| l __ 

(«) While we list the inclination value from HyperLEDA here, our selection of edge-on galaxies was made based on a visual inspection of the 
images. Therefore we have included NGC7814 because of its thin and straight dust lane, but excluded NGC253 because of its visible spiral 
structure. 


Lower helium-burning branches (lowHeB, green): contains 
intermediate-mass core helium-burning stars with ages 
mainly between 100 Myr and 400 Myr, with only little 
contamination by different ages. 

Asymptotic giant branch (AGB, yellow): contains 

intermediate-mass shell helium-burning stars with ages 
mainly between 0.5 Gyr and 2 Gyr, but it might also contain 
a few stars as young as 350 Myr or as old as 6.5 Gyr 
Red giant branch (RGB, red): contains low-mass shell 
hydrogen-burning stars that are mainly older than 3.0 Gyr 
but also younger AGB stars, and even lower HeB stars may 
contribute here. 

Throughout this paper, we use the median ages of these pop¬ 
ulations (see Tablefor all plots showing the ages of the popu¬ 
lations, and their error bars encompass the range between the 10 
percentile and the 90 percentile of the age distributions. 


Completeness correction Not all stars that are in our field can 
be reliably measured. Stars near the detection limit might be 
missed because of statistical variations in the photon flux, and 
stars in crowded or bright sky regions might also be missed. We 
quantified the completeness of our observations through the re¬ 
covery fraction from the artificial star tests. For each popula¬ 
tion the fraction of recovered artificial stars was calculated in 
each spatial bin, that is, we created a completeness map for each 
popul ation. This is different from the approach in the GHOSTS 
DRl ( [Radburn-Smith et al.|2011| ), where the completeness was 
presented as a function of magnitude, color, and sky brightness, 
but not as a function of position. The change in the approach 
was necessary because we found that the completeness can vary 
within a field independent of sky brightness. In particular, we 
found that the completeness depends on the y coordinate of the 
HST images, with a decreasing completeness toward the chip 
gap. This is probably caused by GTE effects, which affect the 
central regions of an image more than the edges. 


Table 2. Properties of the age distribution of the defined stellar 
populations, derived from models assuming a constant SFR and 
a flat metallicity distribution. 



MS 

upHeB 

lowHeB 

AGB 

RGB 

mean(age) [Gyr] 

0.028 

0.089 

0.21 

1.36 

5.07 

std(age) [Gyr] 

0.041 

0.066 

0.19 

0.94 

3.96 

median(age) [Gyr] 

0.010 

0.082 

0.18 

1.16 

4.41 

10 percentile 

0.005 

0.037 

0.12 

0.52 

0.81 

90 percentile 

0.037 

0.146 

0.27 

2.52 

11.5 


2.3. Creating stellar density maps 

We used the star catalogs of the GHOSTS survey (see Sect. |2.1| ) 
to create star c ount density maps for each of the populations 
defined in Sect. |2.2| The raw star count maps are simple two- 
dimensional histograms of the stars in each CMD selection box. 
The bin size of the maps is (7.2")^. These raw star count maps 
were then corrected for incompleteness effects, using complete¬ 
ness maps created from the artificial star tests. They were also 
corrected for masked regions in the original images. We used 
the average of all fields in overlapping regions of different fields. 
Finally, bins with a completeness lower than 0.5 and bins whose 
areas were masked to more than 40% in the original image were 
excluded?. 


2.4. Fitting methods 


To quantify the structural parameters of the disks, we fit dif¬ 
ferent galaxy models to the star count maps. Because the star 
count maps often contain very low numbers of stars per pixel, 
the quality of the fit has to be calculated with a Poissonian like¬ 
lihood estimator instead of the often used; ^^. Such a Po issonian 
likelihood estimator was fir st proposed by Cash ( 1919) and has 
strongly been argued for by jDolphm ( [2002 ). It can be derived in 
the same way as the;^^ statistics, but starting from the Poissonian 
probability function 


n; 

rUi 

ml 


( 1 ) 


with rrii denoting the expected counts in pixel i in the model and 
Hi the observed counts in that pixel. Then the maximum likeli¬ 
hood model can be obtained by minimizing the Poissonian like¬ 
lihood ratio (PLRf] 


PLR = —2 In P = —2 ^ (n, In m, — m, + n,(l — In w,)). (2) 


^ The last term in the sum only depends on the data, which of course 
do not change during the fit. Therefore it can be omitted without effect¬ 
ing the best-fit model, but including this term gives the advantage that 
each term under the sum is greater than zero (which is important for 
some numerical minimization routines) and that for large rif the PLR 
converges to the same value as 
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Fig. 1. Left panel. An artificial CMD generated with a constant star formation rate and a fiat metallicity distribution function (-2.2 < 
[Z] < 0.2). The colored boxes are the population boxes as described in the text. Stars are colored according to their age (see color 
bar). Right panel The age distribution of stars in the five population boxes in the left figure. 


Determining scale heights: To analyze the vertical structure of 
the population at different radii from the galaxy center, we fit an 
isothermal sheet model ( |van der Kruit & Searle||198Ta| to the 
stellar surface density profiles. The model reads 


n{z) = n{) sech^ 


I z-Zc 
\ Zo 




( 3 ) 


where no is the star count density in the midplane, Zc the center 
of the stellar distribution relative to the galaxy midplane (de¬ 
fined through the center coordinates and the position angle of 
the galaxy), zo the scale height, and ntg the surface density of 
contaminants. 

We used an isothermal model even though many observa¬ 
tions found a more peaked function to better fit the vertical light 
profile of galaxies. The isothermal model was preferable because 
we fit our model to distinct stellar populations of well-defined 
ages, for which the isothermal assumption is more justified than 
for the overall stellar content of a galaxy. 


2D fits: To measure scale lengths, scale heights, and break radii, 
we fit a simple model of an edge-on disk with a broken exponen¬ 
tial as radial profile to the 2D maps of star count density: 


n(x, z) = njjg + TZosech^ 


I z-Zc 
\ Zo 


n{x) 


( 4 ) 


with 


/ \X-Xc\ \ 

( \x-xA \ 

^l(rb/hr,o) hr,o ) 


for \x - Xc\ < rij 
for |x - Xcl > r}j ’ 


( 5 ) 


l(x) 


, that i^ 

= Jexp(- 




dy 


OO 


r exp(- 




") lx 

-dr = 2xKi - 
h 

( 6 ) 


This model assumes a constant scale height along the disk. 
We show in Sect. that this assumption is justified. We also 
note that the model in Eq.j^is a single-disk model, meaning that 
it does not contain any thick-disk, bulge, or halo components. 
This will of course limit its applicability; but for the low-mass 
galaxies in our sample, this model is sufficient to describe the 
distribution of stars in each population. 

For the actual 2D fits we used the program IMFIT]^ (Erwin 


2014), which permits using the Poisson likelihood ratio statistics 
for the minimization process. IMFIT is also designed to make 
the inclusion of additional image functions simple, and we have 
extended it with the broken edge-on disk model of Eq. 

The error bars shown in this paper were calculated with 
IMFIT using a bootstrap analysis of the given data. This means 
that all error bars only include statistical uncertainties. 


2.5. Spitzerdata 

For a comparison of the star count maps from GHOSTS with 
integrated light observ ations we used da ta from the InfraRed 
Array Camera (IRAC; [Fazio et ^|2004| ) of the Spitzer Space 
Telescope. The data were red uced within the Spitzer Edge-On 
Disk Galaxies Survey project ( Holwerda et al.|[2006) ). We here 
briefly describe the data and the reduction process. 

The data contain mosaics of 32 edge-on disk galaxies in all 
four IRAC channels. They were taken in a dedicated GO pro¬ 
gram (GO 20268: The Formation of Dust Lanes in Nearby Edge- 
on Disk Galaxies; PI: R. S. de Jong) and were complemented 
with archival data. The observing strategy and data reduction 


where is the x-coordinate of the galaxy center, hrj the inner 
and hr^o the outer scale length, and the break radius. Kfx) is 
the modified Bessel function of the second kind and xKi(x) is 
the projected surface density of an exponential disk seen edge- 


^ This derivation assumes an infinite exponential disk, while we ac¬ 
tually modeled a truncated disk. Therefore the model is not fully self- 
consistent. 

^ http ://WWW. mpe.mpg.de/~ erwin/code/imfit 
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were set up in such a way that the final data products were equiv¬ 
alent in quality to the data products from the Spitzer Infrared 
Nearby Galaxy Survey (SINGS ; [Kennicutt et al.pOOS] ). 

The IRAC Basic Calibrated Data (BCD), together with the 
corresponding uncertainty (BUNG) and indi vidual bad pixel 
mask (BDM SK), were retrieved with Leopard ( [Spitzer Software 
[Team [2008] ). The data reduction was performed with the stan¬ 
dard Spitzer Science Center pipeline for raw I RAC data; the mo¬ 
saics were created with the MQPEX software ([Spitzer Post-BCD 


TooTs Team and Science User Support Team|2012[ ). The reduc- 

tion process includes bias and flatfield corrections, conversion 
from engineering to scientific units, addition of WCS coordi¬ 
nates, combination of single exposures into a mosaic, and it re¬ 
moves cosmic rays. The final images are aligned with the galaxy 
major axes and have a pixel size of 0.75 arcsec. In addition to the 
images, noise maps and mask files were created. 

The IRAC channels 1 and 2, with effective wavelengths of 
X ^ 3.6pm and A ^ 4.5pm, respectively, mainly trace stellar 
emission and have the advantage of being nearly unaffected by 
dust. We used channel 1 to determine the structural parameters of 
the overall stellar population in our galaxies by fitting the edge- 
on disk model (Eq.|^ to the [3.6] data. 

The IRAC channels 3 and 4, with effective wavelengths of 
A ^ 5.1 pm and A ^ 8.0yum, respectively, trace the emission 
of hot dust and polycyclic aromatic hydrocarbon molecules, but 
also include contributions from stellar emission. We used the 
data from channels 1 and 2 to estimate the stellar contribu tion 
to the light in channel 4 (as described in Pahre et al.|2004[ ) and 
subtracted it from the channel 4 data to create images of the non- 
stellar emission in the galaxies. We fit the same edge-on disk 
model as in Eq.|^to these nonstellar images to derive estimates 
of the scale height and -length of the dust in our galaxies. 


3. Results 

In this section we first study the vertical stellar density profiles 
at different galactocentric distances, followed by the radial dis¬ 
tributions at different heights above the midplane. The section 
concludes with a discussion of the observations and a compari¬ 
son with integrated light observations. 


attributed to the extinc ting effects of dust. We discuss this matter 
in more detail in Sect. 14.11 

The radially averaged scale heights that are determined 
through two-dimensional fits are shown as as a function of age 
in Eig. In general, an increase of scaleheights with age is 
visible: While the three youngest populations (MS, upper HeB, 
and lower HeB) have approximately the same scale height, this 
scale height is smaller than the AGB scale height, which is even 
smaller than the scale height of the RGB stars. The relative 
amount of this increase differs from galaxy to galaxy; in IC 5052 
the scale height increases by 100% (from young to AGB) and 
50% (from AGB to RGB), in NGC 5023 by 50% and 33%, and 
in NGC 4244 by only 20% from young to AGB and from AGB 
to RGB. 

We measured the scale heights at different radial positions 
within the galaxies, as can also be seen in Eigs.[^and[^ Eigurej^ 
nicely shows that for most populations the scale heights change 
very little with radial position within each galaxy. While for the 
three young populations the scale heightx are constant along the 
disk, the older populations show some mild flaring. We quanti¬ 
fied the strength of the flaring by fitting a straight line to the scale 
heights as a function of projected radius. The slopes of these re¬ 
gression lines are given in Table 


Table 3. Increase of scale heights with projected radius. 


galaxy 

population 

slope 

[pc/kpc] 

rel. flaring 
[%] 

^max 

[kpc] 

IC5052 

RGB 

22+5 

28+7 

7.1 


AGB 

10+4 

18+7 

7.1 

NGC 4244 

RGB 

7+3 

14+5 

11.9 


AGB 

13+3 

31+9 

11.9 

NGC 5023 

RGB 

9+7 

14+13 

6.8 


AGB 

12+2 

28+6 

6.8 


Notes. Amount of flaring in the older populations of the three galaxies 
(younger populations do not flare). The column ’’slope” gives the ab¬ 
solute flaring dzoldR, the column ”rel. flaring” the relative increase of 
scale height (zo{Rmax)-Zo(R=0))/zo(R=0) over the full observed range in 
radii, and the column R^nax the maximum galactocentric radius to which 
a scale height was measured. 


3.1. Vertical profiles 


We split the observed stars into the five age groups described in 
Sect. [2.2[ The stellar density maps for all of these populations 
are shown in Eig.|^ It is obvious that older populations are more 
extended than young ones. While this is some extent also visible 
in the radial direction, it is clearly evident in the vertical direc¬ 
tion. 

The thickness of the disks can be well examined in Eig. 
Vertical star count profiles are shown there at different radial po¬ 
sitions within the galaxies. The RGB stars (show in red) are the 
dominant population in all our galaxies at all radii; they have 
the highest star count density and the largest vertical extent. To 
quantify the thickness of each population, we fit a sech^ model 
to each radial bin (see Eq.|^; the fits are shown as a dotted line 
in Eig. Remarkably, each profile can be fit well by a single 
sech^ profile (plus background). There is no need to add a sec¬ 
ond component, that is, a thick disk (for more details on this, 
see Sect. 3.3.2[ ). It is also noteworthy that most profiles show a 
clear peak in the center and that only the RGB profiles near the 
galaxy centers have a central dip. Such dips are often observed 
in vertical profile studies (e.g., [Seth et al. 2005b| and are usually 


The slopes in the young populations (MS, upHeB, and 
lowHeB) are consistent with zero flaring. The intermediate and 
old populations in all three galaxies have very similar slopes 
of 10 parsec per kiloparsec, except for the RGB in IC5052, 
which has almost 30pc/kpc. 

3.2. Radial profiles and 2D fits 

We also extracted the radial profiles from our data. The profiles 
at different heights from the midplane are plotted in Eig.[^ While 
the profiles do have some irregularities, they can roughly be de¬ 
scribed as broken exponentials, with steeper slopes outside the 
break. 

Eigure shows that the scale lengths of AGB and RGB 
populations are approximately constant with height above the 
plane. Eor younger populations the height dependence could not 
be measured because there are too few younger stars at large 
heights. 

In two of the galaxies, NGC 4244 and NGC 5023, the aver¬ 
age scale lengths (which were determined through the 2D fits) 
decrease with age, meaning that older populations are more cen¬ 
trally concentrated than younger populations. IC 5052 shows the 
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Fig. 2. Star count surface density maps (logarithmic scale) of stellar populations in three low-mass edge-on galaxies: IC5052, NGC 4244, and NGC 5023 (left to right). Gray 
areas show the regions that were masked because of our crowding limit. Vertical and horizontal gray lines are the bin edges for the extraction of the vertical and horizontal profiles 
in Figs.j^an d0 
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ig. 3. Vertical density profiles of the five populations (red - RGB; yellow - AGB, green - lower HeB, cyan - upper HeB, blue - MS) at different radii in three low-mass edge-on 
alaxies: IC 5052, NGC 4244, and NGC 5023 (left to right). Solid lines are the data, dashed lines are the best sech^ fits. The boundaries of the regions from which the profiles are 
xtracted are shown in Fig.j^as vertical gray lines. 
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Fig. 4. Change of the average scale height with stellar age. The 
blue dashed line shows the results of fits of IC 5052 including 
an additional spheroidal component. The dotted lines are power- 
law zq oc iP fits to the data, excluding the youngest population; 
the power-law indices p are listed on the right. 


opposite trend, and the older populations have fiatter profiles 
than the young populations (see Figs.[7]and[^and Table|^. From 
combining the results for scale length and height, it is clear that 
younger populations are in general thin and radially extended, 
while older populations are thicker and radially more compact, 
as we show in Fig.[^ This is what is seen in the MW for mono¬ 
abundance populations (e.g., Bovy et al.|2012b| ). 



Breaks in radial profiles: The break radii for all populations 
within a galaxy are approximately the same (within a kpc, see 
Fig.[T^, and there is no trend with age. 

Outside the break, the younger populations have steeper pro¬ 
files than older populations. This leads directly to a stronger 
break for younger populations; their profiles are rather fiat inside 
the break and very steep outside, while for older populations the 
scale length changes only little (by about a factor 2-3) across the 
break (see Fig. [n). 


Table 4. Scale heights and scale lengths (inside the break) from 
the 2D fits of the different populations. 


population 

IC5052 

NGC 4244 

NGC 5023 


hr 

zo 

hr 

zo 

hr 

zo 

MS 

1.55 

0.23 

75.36 

0.43 

8.46 

0.21 

upHeB 

1.99 

0.19 

5.63 

0.42 

23.42 

0.20 

lowHeB 

2.50 

0.21 

6.50 

0.44 

19.15 

0.20 

AGB 

2.89 

0.41 

2.29 

0.49 

1.66 

0.31 

RGB 

3.86 

0.61 

2.19 

0.60 

1.84 

0.42 

halo RGB 

- 

2.30 

- 

3.04 

- 

- 

[3.6jj.m] 

2.14 

0.43 

1.84 

0.46 

1.24 

0.30 

PAH [8A/m] 

1.11 

0.38 

1.77 

0.57 

1.08 

0.35 


Notes. All values in kpc. 



Fig. 5. Disk scale heights at different radial position in the three 
low-mass edge-on galaxies IC 5052, NGC 4244, and NGC 5023 
(left to right). Colors are the same as in Fig.|^ 

3.3. Discussion of observed profiles 
3.3.1. Irregularities in the profiles 

IC 5052 This galaxy appears to be lopsided. In Fig.j^the centers 
of all disk fits for IC 5052 are removed by about 1-2 kpc from the 
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Fig. 6. Radial density profiles of the five populations (red - RGB; yellow - AGB, green - lower HeB; cyan - upper HeB; blue - MS) 
at different heights above the plane in three low-mass edge-on galaxies: IC 5052, NGC 4244, and NGC 5023 (top to bottom). Solid 
lines are the data, dashed lines are the best fits. Vertical dotted lines show the break radii of the fits. The boundaries of the regions 
from which the profiles are extracted are shown in Fig. as horizontal gray lines. 
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Fig. 7. Disk scale lengths at different positions above and be¬ 
low the plane in the three low-mass edge-on galaxies IC 5052, 
NGC 4244, and NGC 5023 (left to right). Colors are the same as 
in Fig. 


literature valu^of the galaxy center, which is indicated as the 
zero point of the x-axis in the figure. The AGB and RGB profiles 


^ from HyperLEDA, http: //leda. univ- lyonl. fr/ 


Fig. 8. Change of the average scale length (inside the break) with 
stellar age. The blue dashed line shows the fit results for IC 5052 
including an additional spheroidal component. 


0 1 2 3 4 5 6 0 5 10 15 20 25 0 5 10 15 20 25 



Fig. 9. Scale heights vs. scale lengths (note the different scales 
on the x-axes). 


have their maximum near the literature value of the center, but 
the center of the fit is dominated by the central position between 
the breaks. In Fig.[T^ the contours of youn g stars and RGB star s 
of IC 5052 are shown in an R-band image ( [Meurer et al.|2006) . 
The innermost RGB contours coincide well with the brightest 
region of the R-band image and with the HyperLEDA position of 
IC 5052. In contrast, the center of the outer RGB and the young 
star contours lie farther to the southeast and coincide better with 
the dynamical cen ter of the galaxy, taken from HI observation 
( [Peters et al.|2013| ). 

In the residual images from our two-dimensional fitted mod¬ 
els (see Fig jC.lj ) the overdensity appears approximately circu¬ 
lar. We therefore repeated the fits of the RGB and AGB popula¬ 
tions in IC 5052 with the overdensity described by an additional 
Sersic component and evaluated the structural parameters of the 
disk again. While the scale heights and the break, the radii of the 
two component fits are approximately the same as in the one- 
component fit, the scale lengths are significantly larger. The re- 
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Fig. 10. Change of the break radius with stellar age. The blue 
dashed line shows the fit results for IC 5052 including an addi¬ 
tional spheroidal component. 



Fig. 11. Change of the break strength, i.e., the ratio of inner to 
outer scale length, with stellar age. The blue dashed line shows 
the fir results for IC5052 including an additional spheroidal 
component. 


suits of the two component fits are shown in Figs.|^[^[^ and 
[TT]as a dashed line. 


NGC4244 This galaxy also has an asymmetric radial profile 
(see Fig.[^. At the northeastern (in our plots, at the left) end of 
the disk, all popula tions show a break at about 9 kpc (see also 
|de Jong et al.|2007] ). On the southwestern (right) side the break 
occurs at about 7 kpc, followed by a plateau between 8 kpc and 
12 kpc, where the profile again breaks. These structures cannot 
be modeled by our simple edge-on galaxy models, and thus we 
only used the northeastern part for our fit (up to x < 4 kpc). 

3.3.2. Additional thick-disk or halo component? 

It is now widely accepted that galaxies contain multiple (disk) 
components. In contrast to this, as evident in Fig. the ver- 



Fig. 12. R-band image of IC 5052 (Meurer et al. 


2006). The den¬ 


sity of RGB stars is shown as red contours, the density of young 
stars as blue contours . The yellow reticl e indicates the dynamic 
center of the HI gas ( [Peters et al.|2013| ), the green cross shows 
the central position of the galaxy from HyperLEDA. The black 
boxes show the (approximate) coverage of our GHOSTS fields. 


deal profiles of our low-mass galaxies could be fitted well by 
single-disk models. This is remarkable because more than three 
orders of magnitudes in surface density is covered between the 
central and the outermost parts of the galaxies; this corresponds 
to about eight magnitudes in surface brightness. This is signifi¬ 
cantly deeper than most of the integrated light observations that 
have been used to detect the thick disks. 

However, the profiles presented so far were not optimized 
for low-density regions. As we described in Sect. |2.1| GHOSTS 
has two sets of culls for optimal star selection, one optimized 
for crowded regions, the other optimized for sparse regions. In 
previous sections we used the crowded field culls to optimally 
sample the inner regions of the galactic disks. We now use the 
sparse field culls, which are optimized for low-density regions. 

Additionally, the fine binning of the data we used before 
made it impossible to detect structures below a surface density of 
0.019 arcsec"^ (which equals one star per bin). By using broader 
bins, we can further extend the analysis to fainter regions. In the 
following we use a 200" wide strip along the minor axis of each 
galaxy and adjust the vertcal bin size from 7" at highest den¬ 
sities (near the midplane) to 40" at lowest densities to create 
deeper vertical profiles and search for additional structural com¬ 
ponents. 

The result of this is shown in Fig. Compared to Fig. 
we now reach about one order of magnitude deeper. A faint sec¬ 
ond component can now be detected in IC 5052 and NGC 4244. 
Since we used star count data, this second component cannot be 
an effect of scattered light and extended PSFs, which might lead 
to similar detection in integrated ligh t studies, if not modeled 
carefully ( jde Jong|2008[|Sandin|2014| ). We fit the whole profile 
with two sech^ components. In both galaxies the second compo¬ 
nents are very faint, their central surface star count densities 
are only 0.6% (IC 5052) and 0.08% (NGC 4244) of the main disk 
central surface brightness. As a result of their larger scale heights 
{zQ,ext = 2.30 + 0.69 kpc for IC 5052 and ZQ,ext = 3.04 + 2.29 kpc 
for NGC 4244), the second components dominate at |z| > 2.2 kpc 
and |z| >2.8 kpc, respectively. If we assume that radial structure 
of the extended component is similar to the main disk, we can 
calculate the star count ratio (and thus the stellar mass ratio) of 
the two components. It is M cx which leads to a mass ra¬ 
tio of MextlMthin=0.0\l for IC5052 and for 

NGC 4244. 
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ic5052 



In NGC 5023 we did not detect an extended component, but 
neither did we reach as faint as in the other galaxies. The profile 
of NGC 5023 is based on a single SNAP observation, and its 
CMD reaches only 1.5 mag below the TRGB. This leads to fewer 
detections of stars and thus a higher noise level. Furthermore, the 
vertical extent of our observations is limited to z ~ ±3 kpc (since 
there are no additional fields on the minor axis), which prevents 
detecting an extended component. Thus, we cannot rule out the 
presence of a faint extended component in NGC 5023. 

We note that the choice of sech^ for fitting the second com¬ 
ponent is a mere choice of convenience. It is not meant to imply 
that these components are really disk components. It is possible 
to fit these components equally well with other functions (e.g., 
exponential, Sersic, power law), but the sech^ allows for the most 
direct comparison between the main disk and the second compo¬ 
nent. The possible nature of the second components is discussed 
in Sect.m 




height [kpc] 


Fig. 13. Surface number density profiles of RGB stars in the cen¬ 
tral regions of the three galaxies (shown in red), the dashed lines 
are the best-fit models, the dash-dotted lines show the different 
components, two sech^ functions, and a constant for background 
contamination. The light blue dashed line is the (background- 
subtracted) surface brightness in the Spitzer 3.6/im filter, the 
solid blue line shows the Spitzer data averaged to the same bin 
size as the star count data, and the horizontal gray line plots the 
estimated noise level of the background of the Spitzer data. 


3.3.3. Comparison with Spitzer data 


We here used Spitzer images from channel 1 (at 3.6yum) for 
comparison with the RGB count profiles. For the transformation 
from RGB counts to surface brightness we used the the same ar¬ 
tificial CMDs as in Sect. l2.2l with a constant SFR and a fiat MDF. 
By equating the number of RGB stars in the artificial CMDs with 
the expected flux of the underlying population, we calculated 
the counts-to-fiux conversion factor. This factor can of course 
only be a rough estimate of the true conversion because the ex¬ 
act transformation depends on the star formation history, which 
is unknown for the real galaxies. But as Fig.shows, where the 
Spitzer profiles are shown as light blue and the RGB counts as 
red lines, the transformation matches the profiles quite well. At 
larger heights RGB stars are apparently the dominant contributor 
to the channel 1 light. 

While the Spitzer and RGB profiles match well, there are 
two main differences in a large part of their overlap regions. 

First, the RGB profiles reach much deeper than the Spitzer 
profiles. The Spitzer profiles are limited by a relatively bright 
background. To determine the sky brightness, we measured the 
mean brightness of a large empty region far away from the galac¬ 
tic disk. The corresponding standard deviation in this region was 
used as an esti mate fo r the noise and is shown as a horizontal 
gray line in Fig. 13^^ Star counts are not affected by sky bright¬ 
ness, but they are limited by the number of contaminants, that is, 
by faint foreground stars and unresolved background galaxies. 
The density of these contaminants is so low that our star counts 
can reach an equivalent surface brightness below 30 [3.6yum]- 
mag/arcsec^. 

Second, the Spitzer profiles have a sharper peak in the mid¬ 
plane regions than the star count profiles. The reason for this 
could be that the Spitzer data contain a mixture of all stellar 
populations, as well as dust and PAH emission ( |Meidt et al. 
2012| ). The younger populations, which have a much smaller 
scale height than the RGB stars (see Sect. |3.1| ), will add extra 
light to the central regions of the profile. It could also be that 
we miss some of the RGB stars in the central regions due to 
crowding effects or dust. While we did correct for a decreasing 
completeness with higher star count densities, we did not correct 


The Spitzer images also contain systematic deviations from a con¬ 
stant background level, e.g., large round areas of increased brightness 
in the image of NGC 5023. One of these bright spots lies near the cen¬ 
ter of NGC 5023 and can be seen in the vertical profile as the bump at 
z ^ -2 kpc. 
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for dust effects that will influence our optical data (where the star 
counts are made) more than the infrared Spitzer data. 

In Table 1^ the scale heights and -lengths of the five popula¬ 
tions and the Spitzer data are listed. The Spitzer scale heights lie 
within the range of the scale heights of the resolv ed p opulations 
and are close to the value of the AGB scale heightq^ supporting 
a interpretation as a combination of these different populations. 
The Spitzer scale lengths (of NGC 4244 and NGC 5023), on the 
other hand, are shorter than those of any population, which pre¬ 
cludes an interpretation of the Spitzer data as a simple combina¬ 
tion of the stellar populations. From the two possible reasons for 
this mismatch, incompleteness and/or dust, we asses the former 
as unlikely. Through the extensive artificial star test we have a 
good estimate of the completeness, and we have corrected for 
it. Furthermore, all regions with a completeness below 50% are 
excluded from the analysis, which leads to small uncertainties in 
the final star counts. On the oth er ha nd, we do see signs of dust. 
Its effects are discussed in Sect. 14.11 

3.4. Summarizing results 

Before we move on to the discussion, we summarize all our re¬ 
sults from this section: All three of the analyzed galaxies show 
the following trends in their structural parameters: 

- The scale height of stellar populations increases with age, 
but only on timescales longer than about 300 Myr. 

- The scale height of the young stellar populations is essen¬ 
tially constant along the disk, while AGB and RGB stars ex¬ 
hibit a mild flaring. 

- The scale length of stellar populations is essentially indepen¬ 
dent of the height above the plane. 

- The break radius is independent of stellar population age. 
Two out of three galaxies (NGC 4244 and NGC 5023) have sim¬ 
ilar trends in their radial profiles: 

- The scale length of stellar populations decreases with age. 

- The break strength, that is, the ratio of inner to outer scale 
length, decreases with age. 

A third galaxy (IC 5052) has an increasing scale length with age, 
which leads to a constant break strength. This galaxy also has an 
additional overdensity of old stars, which lies off the center of 
the disk. 


4. Discussion 

4.1. Possible effects of extinction by gas and dust 

In the previous section we have used a sech^ profile for the ver¬ 
tical and a Bessel Ki function for the radial profile of the star 
count surface density distribution of our galaxies. By doing this, 
we have implicitly assumed that the galaxies are transparent, 
meaning that the surface density is the result of the density in¬ 
tegrated along the line of sight through the whole galaxy. The 
presence of dust will lead to a deviation from this idealized as¬ 
sumption. 

Determining the real distribution of dust in external galaxies 
is very difficult and beyond the scope of this paper; so we can¬ 
not correct for extinction effects directly, but we can model the 
effects of different dust distributions on star count maps of an 
idealized galaxy model. 

Table lists the results of a sech^ fit. Since the Spitzer profiles are 
more peaked than the stellar profiles, they are better fit by a generalized 
sech^/^ profile with n > 1. This changes the best-fit Spitzer scale heights 
to be close to the RGB scale height (see Appendix [A|. 


We modeled the density of both the stars and the dust as 
exponential disks with a sech^ distribution in the vertical direc¬ 
tion. We calculated the observed surface brightness by solving 
the standard differential equation for the ID radiative transfer 
along the line of sight. 


dl 

-— = -kI 6, 
dx 


( 7 ) 


with K cx pdusti^. y, z) and e oc p^tarix, y, z). 

It might be questioned whether a radiative transfer equation 
is appropriate for treating star counts. But in fact the effects of 
absorption are similar: It will attenuate and redden the light of a 
star until it falls below the detection limit of our observations (or 
out of our CMD selection boxes); and the higher the absorption, 
the higher the fraction of stars that are removed from the CMD. 
Furthermore, kostenlose email teststar counts have the advantage 
that they are not affected by scattered or reflected light: the light 
of a star that is absorbed will not be re-emitted to become part of 
our star counts; thus the emissivity e is proportional only to the 
local stellar density and has no scattered light term. The details 
of the effects of absorption on star counts do, of course, depend 
on the distribution of stars in the CMD, that is, on the luminosity 
function and, thus, on the star formation history. In this sense, 
our approach is very simplified, but a full extinction modeling is 
beyond the scope of this paper. 

We used three parameters for the dust and stars to model 
the 3D distribution: central density po, scale height zo and scale 
length hr. The density is then given by 


p(x, y, z) = po exp(- + y^lhr)sech^(z/zo). 


( 8 ) 


The parameters of the stellar disk were kept fixed, with hr = I 
and Zo = 0-2. The scales of the dust were varied to be ei¬ 
ther larger or smaller than the stellar ones. We analyzed four 
models with all possible combinations of hr e {0.5,2.} and 
Zo € 10.1,0.4}. The resulting vertical and radial profiles are 
shown in Fig.pl] 

The qualitative effect of the dust on the vertical and radial 
profiles is similar: If the scale height (scale length) of the dust 
is smaller than that of the stars, there is a dip in the vertical (ra¬ 
dial) profile. If the scale height (scale length) of the dust is larger 
than that of the stars, the vertical (radial) profile is flattened in 
the center, but no dip can be seen. The dept h of the dip or the 
flattening depends on the dust mass (see also Seth et al.|2005b 
for a similar analysis)]^ 

We compared the results of our simple dust models (Fig. p^ 
w ith th e observed vertical profiles (Fig. and also Figs. fc.l| 
to |C.3| ). The observed profiles have small dips only in the RGB 
profiles, while the profiles of the AGB and young stars still have 
peaks in the midplane. This suggests that the dust scale height is 
smaller then the RGB scale height, but larger than the AGB scale 
height. This is consistent with our measurements of the scale 
heights of the PAH emission, which is near or above the AGB 
scale height (see Table|^. There is also good agreement between 
our PAH scal e height and the dust sc ale height in NGC 4244 as 
measured by Holwerda et al. ( 2012| ) using far-IR observations 
with Herschel. 

Such large dus t scale heights agree well with the results of 
|Seth et al^ ( |2005b| ), who came to the same conclusion for a sam¬ 
ple of low-mass galaxies (including NGC 4244), and with the 


All four dust models have the same central density. Therefore mod¬ 
els with larger scale heights or scale lengths have higher total dust 
masses. 
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Fig. 14. Effects of dust on ideal galaxy models; Each row shows a model of different dust properties, from top to bottom they are 
defined by (hz 4 ust,hr,dust) = (0.1,0.5) (0.4,0.5), (0.1,2.), and (0.4,2.). Left column: vertical profiles at different radial positions (as 
given in the legend in terms of stellar scale length); middle column: radial profile at different heights above the plane (see legend, 
given in stellar scale heights); right column: face-on view of the position where the optical depth reaches one (r = 1) for different 
heights above the plane (see legend, given in stellar scale heights); the dashed lines are circles of one, two, and three scale lengths 
of the stellar component. 


results of |Dalcanton et al.] ( |2004| ), who found the ISM of low- 
mass galaxies {Vrot < 120km/s) to be more diffuse and vertically 
extended than in massive galaxies. 


the star count data is affected by dust. Thus, the true scale length 
of the underlying populations might be smaller than given here. 


The dust also affects the radial profiles and can lead to a dip 
or a flattening (see Eig. In the data we do not see any central 
dip in our radial profiles (Eig.[^, but a flattening cannot be ruled 
out. Taking into account the smaller scale length of the Spitzer 
3.6/im data, it appears even very likely that the scale length of 


4.2. Direct comparison to other structural measurements in 
these galaxies 


Star count measurements [Seth et alT] ( |2005b| ) analyzed re¬ 
solved stellar populations in the same three galaxies and di- 
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vided them into three age groups (young, intermediate=AGB, 
old=RGB). At that time, they did not have enough data sets to 
study scale lengths or flaring, but they were able to measure scale 
heights. Their results are given in Table they agree with our 
measurements. 


In a series of papers, [Tikhonov & Galazutdinova| studied the 
resolved stellar popula tions of a number of galaxies. [Tikhonov &| 
Galazutdinova| ( [2005| ) analyzed NGC 4244, among other galax¬ 


ies. They used one central ACS fleld (which in our fleld list 
is FieldOl) and two WFPC2 flelds. In agreement with our re¬ 
sults, they found an increasing thickness of the stellar popula¬ 
tions from blue to AGB and RGB stars, but they did not quan¬ 
tify this. In one of the WFPC2 flelds, they claimed to detect 
the transition from a thick disk to a halo in form of a flatten¬ 
ing of the vertical gradient of RGB star density at a height of 
z= 2.7 kpc. This is a simila r heigh t as we found for the transi¬ 
tion to a weak halo in Sect. 3.3.2| but a detailed comparison is 
not possible because their paper lacks a description of the se¬ 
lection of the stars, completeness correction, and treatment of 
unresolved background galaxies. 


[Tikhonov et al.[ ( [2006| ) have analyzed the stellar populations 
of IC 5052 and NGC 5023. They also found an increase of thick¬ 
ness with increasing stellar age and a break in the vertical density 
gradient of RGB stars at z= 1.9 kpc and z= 1.6 kpc for IC 5052 
and NGC 5023, respectively. While for IC 5052 this could be the 
same transition we found in Sect. |3.3.2| we do not see such a 
break in NGC 5023. A more detailed comparison is not possi¬ 
ble due to the only rough description of their data and missing 
flgures in their manuscript. 


Table 5. Scale heights from the literature. 



IC5052 

NGC 4244 

NGC 5023 

Ref. 

counts Zo 





- young 

0.26 

0.33 

0.20 

1 

-AGB 

0.47 

0.44 

0.29 

1 

-RGB 

0.66 

0.55 

0.39 

1 

NIRzo 

0.27/0.47 

0.60 

0.38 

2 


0.25/0.73 



3 


0.19/0.65 



3 



0.60/0.74 


4 

optical Zo 


0.58 

0.46 

5 

dust Zo 


0.6 


6 


Notes. All values in kpc. All scale heights are scaled to our assumed 
distances. 

References. ( 1) |Seth et al!](|2005b| l ;(2)|Salo et al.|(|2013|l; (3)|Comer6n 
let al.|(|2011a|); (4)|Comer6n et al.|(|2011b|); (5)|van der Kruit & Searle 
( |1981a||1982^ ; (6) |Holwerda et al.jpOl^ 


Infrared measurements The galaxies ar e also listed in the S^G 
survey ( [Sheth et~ar 2010 Regan|[2013|), which also includes 
structure decompositions ( Salo et al.[|2013| ). [Salo et ak] ( [2013| ) 
found somewhat larger scale heights (IC5052: 270/740 pc, 
NGC 4244: 600 pc, and NGC 5023: 380 pc) than we did for the 
3.6jum images that range near our RGB scale heights. 

IC 5052 is included in the sample of Comeron et al. ( 2011a ), 
who also used S^G data. They fit a thin (indexed with t) and a 
thick (indexed with T) disk with zo,t = 250 pc and zo,t = 730 pc 
in the outer part and zo,t = 190 pc and zo,t = 650 pc in the in¬ 


ner part of the galax}p^ They obtained a ratio of column mass 
densities = 2.1 in the outer part and = 2.5 in the 

inner part. They stated that their data are also compatible with a 

single-disk structure. _ 

NGC 4244 was separately analyzed in Comeron et al. 
( |2011b[ ). They claimed to have detected a very subtle thick disk 
in NGC 4244, but only on one side of the galaxy, while on the 
other side a single exponential is sufficient to describe the ver¬ 
tical profile. They derived = 350 pc on the northern side and 
hzj = 300 pc and = 370 pc on the southern side, which cor¬ 
responds to zo = 700 pc zo,t = 600 pc and zo,t = 740 pc for a 
sech^ fit. This is significantly larger than our results. A possible 
reason for this discrepancy is the use of different fit functions. 
While an exponential and a sech^ do coinci de for z > 4/z^ = 2zq, 
they differ for lower values. Because both Comeron et al. and 
we also fit regions nearer to the plane, the scale heights for the 
exponential fit are larger than those for the sech^ fit. 


Visual measurement ^van der Kruit & Searle| ( [T981a||1982a| ) 
have measured g-bancp] scale heights as a function of position 
in NGC 4244 and NGC 5023. They found constant scale heights 
along the plane with values of zo = 0.58 kpc for NGC 4244 and 
Zo = 0.46 kpc for NGC 5023. This agrees with our RGB scale 
heights. 


4.3. Disk heating timescales and mechanisms 

The heating of galaxy disks has mostly been studied in the Milky 
Way, mostly due to the lack of good age indicators in exter¬ 
nal galaxies. In the Milky Way the heating is usually examined 
through the age-velocity dispersion relation (AVR), which is 
equivalent to an age-scale height relation, because scale height 
hz and vertical velocity dispersion are connected by 

crl = cnG'Lhz, (9) 


where G is the gravitational constant, X the surface mass density, 
and c a constant that depends on the form of the vertical mass 
profilj^ 

Even though the AVR has been studied for decades, there is 
still no consensus on the form of the AVR. Some authors found a 
smooth power-law increase of the velocity dispersion with time 
( [Wielen|1977t [Nordstrom et al.|2004t[Kovar et al.|2009|), while 

others found a flattening for stars older than a few Gyr ( Carlberg| 


very old stars (> 9 Gyr) in the thick disk ( 

Edvardsson et al. 19931 

Quillen & Garnett|2001l 

Anguiano et al. 

201 21 1 Anguiano|2012). 

Possible reasons for the c 

iscrepancies might be difficulties in age 


ing all stars vs. excluding moving groups and thick-disk stars). 
But the authors’ choice of model can likewise influence their 


conclusions. For example, Seabroke & Gilmore (2007) showed 


The original values in their paper are different (180 pc, 530 pc, 
140pc, and 470pc); we have converted them to our scale: While they 
used exponential s cale heights, we used sech^, which differs by a factor 
of 2. Furthermore, jComeron et al.| assumed a larger distance (8.1 Mpc). 
We have corrected for this. 

The authors denote their band ”J”, but w ith an effective bandpass 
from 400 nm to 540 nm { [van der Kruit] 197"^ it lies between the B and 
V band and covers a wavelength range similar to the SDSS g-band. 

e.g., for a sech^ it is c=2 and for an exponential c=3/2. 

For example, an age uncertainty of only 30% for the oldest stars, 
which is much lower than the uncertainties in current data, can er ase a 
sharp increase in velocity dispersion at 9 Gyr ( jMartig et al.|2014^ . 
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that although the data of Nordstrom et al. ( |2004 ) are indeed con¬ 
sistent with a power-law increase of cr^, neither can they reject a 
saturation of cr^ at ages older than 4.5 Gyr. 

As possible mechanisms for heating the disk, scattering by 
giant molecular clouds (GMC) w ere the first to be proposed by 
[Spitzer & Schwarzschild 


f l951| ), followed by spiral structure 
a nd massive compact ha lo objects 


(MACHOS, e.g., black holes; [Lacey & Qstriker|1985] ). None of 
these alone were found to be supported by the observations (see, 
e.g., [Binney & Tremain^|2008| and reference s therein). Other 
proposed mechanism include minor mergers ( [Toth & Ostriker j 
|1992j ), which can lead to a sudden increase of velocity disper¬ 
sion, when the satellite is relatively massive ( [Quinn et al.|1993| ), 
or to a smooth heating, when many small galaxies are accreted 
over a long time span ([Velazquez & White |1999[ ), resonances 
from the bar ( [Kalnaj s[ 1991 [ ) and dissipating star clusters ( [Kroupa| 
2002[ ). We note that the two most accepted models, namely scat¬ 


tering by GMC and spiral arms, would lead to a saturation. This 
is easily understood for the vertical motion cr^, since GMC and 
spiral arms act on stars in the very midplane of the disk, and 
any star with a high vertical velocity will spend most of its or¬ 
bital time far away from the plane. Thus, it would no longer be 
affected by GMCs or spiral arms. 

Our data do not have the necessary age resolution to answer 
the question of a possible saturation of the heating, but we can 
constrain heating efficiencies. We fit a power law zo ^ ^ to 
our data (see Fig. [^ excluding the data point of the youngest 
stars) and found indices {3 in the range 0.09-0.32. According to 
Eq. [^ this corresponds to a power law for the velocity disper¬ 
sion (Tz t>c with a = pjl. This is significantly lower than the 


Milky Way value 


> 0.3 (e.g., [Hanninen & Flynn[[2002l 

1 but consistent with the values for low- 

Nordstrom et al. 2004[ 

n 

mass galaxies found 1 
nation for the lower di 

)y [Seth et al. ( 

2005b). A probable expla- 

isk heating in 

ow-mass galaxy is that the 


two main heating mechanisms, scattering by GMCs and by spiral 
structure, are less efficient because low-mass galaxy are gravita¬ 
tionally more stable, resulting in a less strong spiral structure, 
and a more extende d, but less dense interst ellar medium, result¬ 
ing in fewer GMCs ( [Dalcanton et al.|2004| ). 


4.4. Structural evolution of stellar populations 


[Martig et~ST ( 2014a|b[ ) have analyzed seven simulations of disk 
galaxies representing a variety of merger histories. They have 
separated each galaxy into so-called mono-age populations, each 
representing an age range of 500 Myr, and analyzed how their 
structural parameters (scale length and -height, velocity disper¬ 
sion) change as a function of age (at z=0). 

In galaxies with a quiescent merger history they found a 
smooth decrease of scale length with age and a smooth increase 
of scale height (and vertical velocity dispersion). Only for very 
old stars (age> 8.5 Gyr), a jump in vertical velocity dispersion 
can be seen. These very old stars are centrally more concen¬ 
trated and show strong flaring; they are more likely associated 
with a bulge or inner halo component. The physical reason for 
the smooth evolution of scale lengths is seen in radial migration. 
Stars are born in an irregular, but relatively fiat profile with a 
break. With time, stars migrate and populate the outer disk, thus 
erasing the break. At the same time, the central parts contract 
(due to conservation of angular momentum), which leads to a 


steepening of the radial profile. The vertical heating in their sim¬ 
ulation is caused by spiral arms, overdensities in the disk, the bar 
and bending waves, but not by radial migration. 

As expected, mergers have a significant effect on the evolu¬ 
tion of scale height: A merger leads to a sharp increase in the 
AVR by increasing the velocity dispersion or scale heights of all 
populations born before the merger. But it also has an interesting 
effect on the scale length evolution: Near the time of a merger, 
the normal trend of scale lengths with age is reversed. For a cer¬ 
tain age range (usually about 2 Gyr), the scale length increases 
with age. This leads to V-shaped structures in a scale length vs. 
age plot. During this period, scale height and scale length are 
correlated. 

The increase (decrease) of scale height (scale length) with 
age and the corresponding anti-correlation in a h^-hr plot that 
Martig et al.[see in quiescent galaxies is very similar to what we 
see in NGC 4244 and NGC 5023 (see Fig.^. 


The exceptional behav ior of IC 505 2 can also be explained in 
the light of the results of Martig et if we assume that IC 5052 
recently experienced a merger. This assumption is strongly sup¬ 
ported by the irregular radial profile of IC 5052 (see Sect. [3.3.1] ); 
the overdensity of old stars that is found with an offset from the 
dynamical center of the galaxy could be the remnant of a merg¬ 
ing satellite. Furthermore, the residual map from the fit (Fig. [C.l[ ) 
reveals a stream-like feature that connects the overdensity on the 
northwestern (in our plots left) side of the disk with a smaller 
overdensity at the southeastern (right) side. [Martig et al.| ( [2014a| ) 
have also shown that a merger will lead to a strong flaring of the 
pre-existing disk. That fits the fact well that the old population 
of IC 5052 has the strongest hints on flaring (see Table [^ and 
the highest heating efficiency (see Fig.[^ of all our galaxies. 


Bird et al!] ( |2013| ) have analyzed a high-resolution cosmologi¬ 


cal simulation of a Milky Way-mass galaxy. They also separated 
their simulated galaxy into different age groups and found an 
increasing scale height an d a decreasing scale le ngth with age. 
These trends are similar to Martig et al. ([2Q14a Q and our obser¬ 
vations, but in contrast to Martig et al. ( 2014a|b ), these trends are 
imprinted on the populations during their formation and change 
only little during the evolution of the galaxy. In their model the 
age evolution of the structural parameters is a direct consequence 
of the time evolution of the properties of the gas from which the 
stars are formed. While the gas at high redshift has a large verti¬ 
cal extent and is mainly dispersion supported, it becomes thinner 
and radially more extended with time; and the stars follow this 
behavior. 


Another feature in the galaxy of [Bird et al.| ( [2013} can also 
be compared with our observations. Their galaxy has a strong 
flaring in all age groups and velocity dispersion profiles that are 
almost constant with radius for each age group. This stro ngly 
disagrees with our observations and with integrated light ([van 


der Kruit & Searle|1981a|blll982a|bllShaw & Gilmore|1990llde 

Grijs & van der Kruitjl996[[^^ 

and kinematical observations ([van 

der Kruit & Freeman|| 19841 

19861 |Bottema|| 19931 |Martinsson 


et al.[[2013[ ), whichTdl found constant scale heights and expo¬ 


nentially decreasing veloci ty dispersions. The strong flaring re¬ 
ported by [Bird et al.| ( [2Q13| is not a consequence of the evolution 


There is no clear consensus on the value of a for the MW, but almost 
all authors found va lues in the range 0.3 < cr < 0.6, with the exception 
of Anguiano ( 2012| ), who found a = 0.12 for the MW thin disk. 


Some galaxies in these papers do show a hint of an increase of scale 
height at the very edges of the disks. See also |de Grijs & Peletier| ( |199^ 
for a sample of flaring galaxies. 
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or interactions, but it is also already imprinted on the gas at all 
times. 


The small amount of flaring that we observe in the three 
galaxies also contradicts most simulations that model disk for¬ 
mation in a cosmological context . The ubiquity of merger events 
inevit ably produces flaring disks ( [Bournaud et al.|2009 Qu et al 


|2Q11| ). The absence of (strong) flaring could thus indicate a 
rather quiet merger history of these galax ies. A possible alterna - 
tive explanation was recently proposed by Minchev et al.| ( |2015] ), 
who found that the total stellar density can appear unflared even 
when all the mono-age populations flare. This may occur be¬ 
cause the scale lengths change with age and the flare of each 
population occurs at a different radii. The combined flares of all 
populations can then be described as a nonflaring thick disk with 
a radial age gradient. While we do resolve the stars in popula¬ 
tions of different ages, the wide age range especially of the RGB 
component could smear out the flaring effects; but in this case, 
a separate thick-disk component should be present. The absence 
of such a thick disk renders this scenario unlikely for our obser¬ 
vations. 


[Stinson et al.| ( |2013| ) compared a cosmological simulation of 
a Milky Way-like galaxy including detailed chemical model¬ 
ing with the results of jBovy et al.| ( |2012b| ). They also separated 
their model galaxy into mono-abundance populations and found 
an anticorrelation between scale height and scale length. They 
showed that the mono-abundance populations are good approx¬ 
imations of mono-age populations (usually with Icr age spreads 
of less than 1 Gyr) and that the anticorrelation is indeed caused 
by an increase of scale height and a decrease of scale length 
with age. This therefore fully agrees with our observations of 
NGC4244 and NGC5023. [Stinson et al.j ( [2Q13[ ) have neither 
studied the radial variations of scale heights, nor the physical 
origins of these relations. 


4.5. Thin- and thick-disk dichotomy? 


Studies of the age evolution of galactic disks beyond the MW 
(or even the Local Group) similar to the theoretical work in 
Sect. 4.4 were very difficult until recently, and only a few authors 
have trie d it (e.g., [Seth et al.|2005b} [Tikhonov & Galazutdinova[ 
[2012a|b[ ). Thus, most observational data on disk evolution comes 
form thin- and thick-disk decomposition 

Thick disks a re usually found to be thicker by a factor of at 
least >1.6 ( [Yoachim & Dalcanton|2006|), in most cases it 
is even ztI zt > 2.6 ( [Comeron et al.[[2Ql lat [Pohlen et al.[[20Q4[ 
and references therein). The scale lengths of the thick compo¬ 
nents are, wit h only a few exceptions, also larger than those of 
the thin disk ([Pohlen et al.[[2004t [Yoachim & Dalcanton[[2QQ6| 


Comeron et al.|2QlTa ). On a verage, no flaring was f ound in ei¬ 

ther the thin or the thick disk ( [Comeron et al.|2Ql la[ ). 


The question then is how our disks would fit into these sam¬ 
ples and whether our old disks are comparable to thick disks. 

Figure shows the scale height s compared to thin- and 
thic k-disk scale heights reported by [Comeron et al.[ ( [201 la[ ) 
and [Yoachim & Dalcanton[ ( |2006| ). While all our measured disk 
scale heights are consistent with thin-disk scale heights, the 
scale heights of the faint extended components in IC5053 and 
NGC 4244 are larger than typical thick-disk scale heights. 


or from color gradients, but degeneracies between age, metallicity, 
and extinction effects make an interpretation of these gradients ambigu¬ 
ous. 
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Fig. 15. Comparison of our measu red scale heights with dat a 
from [Yoachim & Dalcanton[ ( [2006[ ) and [Comer6n et al.| ( [2011a[ ). 
The colored circles show our results for the different populations 
(color code as before); squares show the faint extended compo¬ 
nent (Sect. [3.3.^ . 


We found old-to-young scale height ratios of about 3, 1.5, 
and 2 (for IC5052, NGC 4244, and NGC 5023, respectively). 
While for IC5052 and NGC 5023 this agrees well with ear¬ 
lier thick-disk studies, it would make the disk of NGC 4244 the 
thinnest thick disk (relative to the thin disk) seen so far. These 
thickness ratios also agree well with the velocity dispersion mea¬ 
surement within the thin disk of the Milky Way (see Sect. [4.3[ ): 
The velocity dispersion of stars increases from less than 10 km/s 
at an age of 1 Gyr to more than 20 km/s at ages of about 6 Gyr. 
According to Eq.[^ such a doubling of velocity dispersion would 
result in an increase of scale height by a factor of four, which 
is already more than we observe. Furthermore, in two of three 
galaxies (namely NGC 4244 and NGC 5023) the old and thicker 
component has a much shorter scale length than the young pop¬ 
ulations, in contrast to earlier thick-disk observations in external 
galaxies. 

Therefore, we conclude that the old disks observed here are 
really the old thin disks of these galaxies. More extended com¬ 
ponents are much fainter and c annot contain more than 2% of 
the stellar mass (see Sect. [3. 3. 2| ). 


4.6. Breaks 


In each galaxy we found the break radius to be the same for all 
populations. The breaks are sharper, that is, the ratio /tr,//^r,o is 
higher, fo r younger population s. T his agrees with earlier obser - 
vations by [de Jong et al.[ ( [2QQ7] ) and [Radburn-Smith et al. ( [2Q12[ ). 
The conclusion that a constant break radiu s favors a dynami- 
cal origin of the bre ak ( de Jong et al.[[2007[ ) was questioned by 
[Roskar et al.[ ( [200^ . In their simulation, the break was caused 
by a decrease in the gas density and the break radius moved out¬ 
ward with time, while gas with higher angular momentum was 
accreted. Despite this change of break radius with time, at the 
end of the simulation the break radius was the same for popu¬ 
lations of all ages. This is due to the radial migration of stars: 
The oldest stars that have the most time to migrate move farthest 
outward, shifting their break to larger radii and weakening their 
strength. 
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A similar behavior is seen in the test particle simulations of 
[Elmegreen & Struck ( |2Q13| ). They started with a flat radial den¬ 
sity profile of stars with a sharp cutoff at a given radius. This 
initial disk was perturbed by transient mass clumps in the disk, 
which led to a double exponential profile with decreasing break 
strength, until at the end, after several Gyr, the break was erased 
and a single exponential had formed. During this process, the 
break radius also moved outwardP^ 

Our findings are in contrast with those of Bakos et al.| ( |2008) ), 
who argued that the break is only a change of stellar populations, 
but not a change in the mass density profile. While we also see 
a change in stellar populations, namely a change in the ratio of 
young to old stars, which has a maximum near the break, it is 
also clear from the existence of the break in all populations that 
the mass profile is also broken. 


5. Conclusions and summary 

We have mapped the stellar populations of three low mass, edge- 
on disk galaxies: IC 5052, NGC 4244, and NGC 5023. Through 
analyzing their CMDs, we separated the stellar content of the 
galaxies into five populations of distinct ages. We measured 
the structural parameters (scale heights, inner and outer scale 
lengths, break radii, and flaring) for each population and ana¬ 
lyzed their dependence on stellar age. 

Each component of each galaxy was well fit with a single¬ 
disk model with constant scale height and a broken exponential 
radial profile. None of the populations showed a need for second 
(e.g., thick) disk component. Along the minor axes of IC5052 
and NGC 4244 we detected very faint additional components in 
the old population, which contained at most 1 % of the main disks 
mass and are thus much fainter than a typical thick disk. These 
extended components are probably the faint halos of their galax¬ 
ies. 

We estimated the amount of flaring in all the populations. 
While the younger populations do not flare, the intermediate and 
old populations show a small increase of scale height by about 
10 pc per kpc in radius. This adds up to a relative flaring of not 
more than 30% over the whole observed radial range, which 
reaches well beyond the break in the radial profiles. While such 
a low amount of flaring is consistent with most earlier obser¬ 
vations, it challenges cosmological simulations of galaxy disks, 
where interactions and mergers usually lead to a stronger flaring. 

All galaxies have an increase of scale height with age, but 
the disk heating efficiency, expressed through the index a in the 
power law of the form cr^ oc oc , is lower than in the 
Milky Way. This can be well understood through the weaker spi¬ 
ral structure (i.e., higher gravitational stability) and more diffuse 
interstellar medium of low-mass galaxy disks. 

The break in the radial profiles appears at the same radius for 
all ages in each galaxy. In two of the galaxies (NGC 4244 and 
NGC 5023), the break strength decreases with age, which sup¬ 
ports break formation models that induce a sharp cutoff of star 
formation (either through star formation thresholds or limited 
gas supply), which are then softened by scattering and migration 
of stars over time. 

In IC 5052 we detected an additional overdensity of old stars, 
which lies off the center of the disk, and a stream-like feature 
in the fit residuals. These were interpreted as an indication of a 


The authors stated that the break radius of the double exponential 
profile was comparable to the initial cutoff, but a detailed look at their 
Fig. 2 shows an increase of break radius with time. 


recent merger. The effects of this merger can be seen in the in¬ 
creasing scale length with age, the larger amount of flaring, and 
the higher heating efficiency in IC 5052 compared to the other 
two galaxies. 
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Appendix A: Scale heights in the Spitzer [3.6jj] data 


We show here the measurements of scale heights along the major axis of the galaxies for the Spitz er data . We fit both a simple sech^ 
and generalized sech^/". The Spitzer profiles are usually more peaked than a sech^ (see also Sect. 3.3.3 ), and we found n > 1. The 
scale heights of the sech^/" fits agree well with the RGB scale heights, and there is no or very little flaring. 

The parameter n of the sech^/" increases toward the center of the galaxies, meaning that the profiles in the inner disks are more 
peaked than in the outer disks. 




Fig. A.l. Upper panels: Scale heights in Spitzer IRAC channel 1 as a function of position along the major axis (from left to right: 
IC 5052, NGC 4244, NGC 5023). Lower panels: parameter n of the generalized sech^/" fits. 


Appendix B: Color-magnitude diagrams 

We present here a number of color magnitude diagrams to show the population selection boxes and to illustrate the variety of CMD 
features, stellar densities, and observation depths. For each galaxy we show one CMD from the central disk, the outer disk, and the 
halo (if available). 
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Fig.B.l. Color-magnitude diagrams from three fields in IC 5052. We show a field from the central disk (left), the outer disk (middle), 
and the halo (right). Dense regions in the CMD are plotted as a Hess diagram with con tour s at 100, 200, 400, and 800 stars per 


(0.1 mag)^. Colored boxes are the selection boxes for the populations defined in Sect, 
completeness level of the observations. 
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and the red dashed line is the 50% 


Appendix C: 2D Fits 

We present here the details of the two-dimensional fits to the star count maps. For each galaxy we show for each population the data 
and the model image. We also show the weighted residuals. We defined the residual wdevi of each pixel in the context of Poisson 
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Fig.B.2. Color-magnitude diagrams from three fields in NGC 4244: Details as in Fig 


B.l 
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(for NGC 5023 no halo field is available). 


Fig.B.3. Color-magnitude diagrams from three fields in NGC 5023: Details as in Fig 


B.l 
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noise dominated data to be 

wdevi = sign(7Z/ - m/) ^J-2(ni Inm/ - m/ -l- 7Z/(1 - In rit)), (C.l) 

with Hi the pixel value of the data and m/ the pixel value of the model. This is analog to the normal weighted residual (rii - mi)lcri 
in the sense that it is positive for rif > m/ and negative for rit < and it is PLR = 
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Fig. C.l. Top row: 2D fit maps of the five populations in IC5052 (from left to rig ht: RGB, AGB, lowHeB, upHeB, MS). Each panel contains three maps: top: stellar density 
maps, middle: best-fit model, bottom: weighted residual (data-model, see Eq. |C.l) . The gray lines show the regions from which the minor and major axis profiles are extracted. 
The two overdensities in the residual map of the fit to the RGB data (left panel); bottom row (middle and right): major and minor axis profiles of the data (solid lines) and the 2D 
models (dashed lines); dotted vertical lines show the break radii. The model lines may show irregularities (as the real data do) due to masked-out regions. 
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Fig. C.2. Same as Fig.|C.lj but for NGC 4244. We restricted the fits to the left side of the galaxy because of irregularities at the right end. 
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